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Abstract 

The ellipsoidally symmetric extension of Buda-Lund hydrodynamic 
model is shown here to yield a natural description of the pseudorapidity 
dependence of the elliptic flow V2(n), as determined recently by the PHO- 
BOS experiment for Au+Au collisions at ^/sjvjv = 130 and 200 GeV. 
With the same set of parameters, the Buda-Lund model describes also 
the transverse momentum dependence of vi of identified particles at mid- 
rapidity. The results confirm the indication for quark deconfmement in 
Au+Au collisions at RHIC, obtained from a successful Buda-Lund hydro 
model fit to the single particle spectra and two-particle correlation data, 
as measured by the BRAHMS, PHOBOS, PHENIX and STAR collabo- 
rations. 



"Nuclei, as heavy as bulls, 
through collision, generate 
new states of matter. " 

(T. D. Lee) 



1 Introduction 

Ultra-relativistic collisions of almost fully ionized Au atoms are observed in four 
major experiments at the RHIC accelerator at the highest currently available 
colliding energies of ^/snn = 200 GeV to create new forms of matter that existed 
before in Nature only a few microseconds after the Big Bang, the creation of 
our Universe. At lower bombarding energies at CERN SPS, collisions of Pb 
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nuclei were studied in the y / ijv r /v = 17 GeV energy domain, with a similar 
motivation. If experiments arc performed near to the production threshold of a 
new state of matter, perhaps only the most violent and most central collisions are 
sufficient to generate a transition to a new state. However, if the energy is well 
above the production threshold, new states of matter may appear already in the 
mid-central or even more peripheral collisions. Hence the deviation from axial 
symmetry of the observed single particle spectra and two-particle correlation 
functions can be utilized to characterize the properties of such new states. 

The PHENIX, PHOBOS and STAR experiments at RHIC produced a wealth 
of information on the asymmetry of the particle spectra with respect to the re- 
action plane ^ El EH EH 03 EH , characterized by the second harmonic moment 
of the transverse momentum distribution, frequently referred to as the "elliptic 
flow" and denoted by v 2 . This quantity is determined, for various centrality se- 
lections, as a function of the transverse mass and particle type at mid-rapidity 
as well as a function of the pseudo-rapidity 77 = 0.5 log( ctz^T )• Pseudorapid- 
ity measures the zenith angle distribution in momentum space, but for particles 
with high momentum, |p| » E\ p \, it approximates the rapidity y = 0.51og( ^^ ) 
that characterizes the longitudinal momentum distribution and transforms ad- 
ditively for longitudinal boosts, hence the rapidity density dn/dy is invariant 
for longitudinal boosts. The PHOBOS collaboration found 0, that V2(r)) is 
a strongly decreasing function of which implies that the concept of boost- 
invariance, suggested by Bjorken in ref. [7] to characterize the physics of very 
high energy heavy ion collisions, cannot be applied to characterize the hadronic 
final state of Au+Au collisions at RHIC. A similar conclusion can be drawn from 
the measurement of the inhomogeneous (pseudo)rapidity dn/drj and dn/dy dis- 
tributions of charged particle production at RHIC by both the BRAHMS 8 
and PHOBOS 9 collaborations, proving the lack of boost-invariance in these 
reactions, as dn/dy ^ const at RHIC. Although many models describe success- 
fully the transverse momentum dependence of the elliptic flow at mid-rapidity, 
v-zipt, T] = 0), see ref. [G3] for a recent review on this topic, to our best knowledge 
and an up-to-date scanning of the available high energy and nuclear physics lit- 
erature, no model has yet been able to reproduce the measured pseudo-rapidity 
dependence of the elliptic flow at RHIC. 

Hence we present here the first successful attempt to describe the pseudo- 
rapidity dependence of the elliptic flow v 2 (rj) at RHIC. Our tool is the Buda- 
Lund hydrodynamic model |12l which we extend here from axial to ellip- 
soidal symmetry. The Buda-Lund hydro model takes into account the finite 
longitudinal extension of the particle emitting source, and we show here how 
the finite longitudinal size of the source leads naturally to a 1*2 that decreases 
with increasing values of \rj\, in agreement with the data. We describe simul- 
taneously the pseudorapidity and the transverse momentum dependence of the 
elliptic flow, with a parameter set, that reproduces the single-particle trans- 
verse momentum and pseudo-rapidity distributions as well as the radius param- 
eters of the two-particle Bose- Einstein correlation functions, or HBT radii, in 
case when the orientation of the event plane is averaged over. All these benefits 
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are achieved with the help of transparent and simple analytic formulas, that are 
natural extensions of our earlier results to the case of ellipsoidal symmetry. 



2 Buda-Lund hydro for ellipsoidal expansions 

The Buda-Lund model is defined with the help of its emission function S(x,p), 
where x = (t, r x ,r y ,r z ) is a point in space-time andp = (E,p x ,p yi p z ) stands for 
the four-momentum. To take into account the effects of long-lived resonances, 
we utilize the core-halo model |14j . and characterize the system with a hydro- 
dynamically evolving core and a halo of the decay products of the long-lived 
resonances. Within the core-halo picture, the measured intercept parameter 
A* of the two-particle Bose-Einstein correlation function is related \TQ to the 
strength of the relative contribution of the core to the total particle production 
at a given four-momentum, 

S(x,p) = S c (x,p) + S h (x,p), and (1) 
S c (x,p) = y/X»S(x,p). (2) 

Based on the success of the Buda-Lund hydro model to describe Au + Au colli- 
sions at RHIC [151 Pb + Pb collisions at CERN SPS (THj and h + p reactions 
at CERN SPS (171 118) . we assume that the core evolves in a hydrodynamical 
manner, 

S c {x,p)d x = —-— — , 3) 

(2TT) i B(x,p) + s q 

where g is the degeneracy factor (g = 1 for identified pseudoscalar mesons, g = 2 
for identified spin=l/2 baryons), and p fJ, d 4 Y l/1 (x) is a generalized Cooper-Frye 
term, describing the flux of particles through a distribution of layers of freeze- 
out hypersurfaces, B(x,p) is the (inverse) Boltzmann phase-space distribution, 
and the term s q is determined by quantum statistics, s q = 0, —1, and +1 for 
Boltzmann, Bose-Einstein and Fermi-Dirac distributions, respectively. 

For a hydrodynamically expanding system, the (inverse) Boltzmann phase- 
space distribution is 

W - 2$) ■ < 4 » 

We will utilize some ansatz for the shape of the flow four- velocity, u v (x), chem- 
ical potential, fi(x), and temperature, T{x) distributions. Their form is deter- 
mined with the help of recently found exact solutions of hydrodynamics, both in 
the relativistic [12112011111 and in the non-relativistic cases |22l 1231 124| , with the 
conditions that these distributions are characterized by mean values and vari- 
ances, and that they lead to (simple) analytic formulas when evaluating particle 
spectra and two-particle correlations. 

The generalized Cooper-Frye prefactor is determined from the assumption 
that the freeze-out happens, with probability iJ(r)dr, at a hypersurface char- 
acterized by t = const and that the proper-time measures the time elapsed 
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in a fluid element that moves together with the fluid, dr = u fl (x)dx fl . We 
parameterize this hypersurface with the coordinates (r x ,r y ,r z ) and find that 
d 3 S^(x|r) = u fJ, (x)d 3 x/u°(x). Using d t r\ r = u°(x) we find that in this case the 
generalized Cooper-Frye prefactor is 

p"d%(a:) = p^u^{x)H(r)d 4 x, (5) 

This finding generalizes the result of ref. [2SJ from the case of a spherically 
symmetric Hubble flow to anisotropic, direction dependent Hubble flow distri- 
butions. 

From the analysis of CERN SPS and RHIC data [H3 USJ gSj , we find that 
the proper-time distribution in heavy ion collisions is rather narrow, and H(t) 
can be well approximated with a Gaussian representation of the Dirac-delta 
distribution, 

H W - (2^2)1/2 (~ I 1a3 L ) ' (6) 

with At -C tq. 

Based on the success of the Buda-Lund hydro model to describe the axially 
symmetric collisions, we develop an ellipsoidally symmetric extension of the 
Buda-Lund model, that goes back to the successful axially symmetric case |12l 
1131 1161 1151 |2U if axial symmetry is restored, corresponding to the X = Y and 
X = Y limit. 

We specify a fully scale invariant, relativistic form, which reproduces known 
non-relativistic hydrodynamic solutions too, in the limit when the expansion 
is non-relativistic. Both in the relativistic and the non-relativistic cases, the 
ellipsoidally symmetric, self-similarly expanding hydrodynamical solutions can 
be formulated in a simple manner, using a scaling variable s and a corresponding 
four- velocity distribution u M , that satisfy 



d^s = 0, (7) 

which means that s is a good scaling variable if its co-moving derivative van- 
ishes ^HEO]- This equation couples the scaling variable s and the flow velocity 
distribution. It is convenient to introduce the dimensionless, generalized space- 
time rapidity variables {f] x ,r]y,r] z ), defined by the identification of 

X Y Z 

(sinhrfc, sinh%, sinhrfe) = [t x — , r y — , r z — ). (8) 

Here (X, Y, Z) are the characteristic sizes (for example, the lengths of the ma- 
jor axis) of the expanding ellipsoid, that depend on proper-time r and their 
derivatives with respect to proper-time are denoted by (X,Y,Z). The distri- 
butions will be given in this r/i variables, but the integral-form is the stan- 
dard d 4 x — dtdr x dr y dr z , so we have to take a Jacobi-determinant into account. 
Eq. (JJJ is satisfied by the choice of 

cosh r] x — 1 cosh rj y — 1 cosh rj z — 1 



vl v2 72 



(9) 



(7, sinhffc, sinh?7 a , sinh?7 z ), (10) 
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and from here on (Xf,Yf,Zf) = (X(tq), Y{tq), Z(tq)) — (Xi 7 X 2 ,X 3 ), assum- 
ing that the rate of expansion is constant in the narrow proper-time interval of 
the freeze-out process. The above form has the desired non-relativistic limit, 



r 2 r 2 
' x . y 



2Xj 2Yf 2Zf 



(11) 



where again (X f ,Y f ,Z f ) = (X(r ), Y{tq), Z(t )) = (X 1 ,X 2 ,X 3 ). From now 
on, we drop subscript /. From the normalization condition of u /i (x)M At (a;) = 1 
wc obtain that: 



7 = y 1 + sinh r\ x + sinh r\ y + sinh rj Zl (12) 

For the fugacity distribution we assume a shape, that leads to Gaussian profile 
in the non-relativistic limit, 

T(x) T 

corresponding to the solution discussed in refs. [22123 EH- We assume that 
the temperature may depend on the position as well as on proper-time. We 
characterize the inverse temperature distribution similarly to the shape used in 
the axially symmetric model of ref. ^1 EH > an d discussed in the exact hydro 
solutions of refs \'22\ 123) . 

1 ^l + ^.)fl + ^fe#) (14, 



T(x) T V T s J V T e 2A 

where To, T s and T e are the temperatures of the center, and the surface at 
the mean freeze-out time To, while T e corresponds to the temperature of the 
center after most of the particle emission is over (cooling due to evaporation 
and expansion). Sudden emission corresponds to T e = To, and the At — > 
limit. It is convenient to introduce the following quantities: 



2 = Tp-T s = I AT\ 
Ts \T A 

2 _ T — T e _ / AT\ 
Te " \ r A 



(15) 



(16) 



In the above approach we assume the validity of the concept of local thermal- 
ization and the concept of ellipsoidal symmetry at the time of particle produc- 
tion. We do not know exactly, what freeze-out condition is realized in Nature. 
Our above formulas can be also considered as a general, second order Taylor ex- 
pansion of the inverse temperature and logarithmic fugacity distributions, while 
maintaining ellipsoidal symmetry. We attempt to determine the coefficients of 
this Taylor expansion from the data in the subsequent parts. As the saddle-point 
calculation presented below is sensitive only to second order Taylor coefficients, 
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any model that has similar second order expansion leads to similar results. A 
more theoretical approach is to solve relativistic hydrodynamics for ellipsoidally 
expanding fireballs and to apply the presently most advanced freeze-out criteria, 
for example the method of escaping probabilities developed by Akkelin, Hama 
and Sinyukov in ref. 

We also note that the applied distribution of freeze-out temperatures goes 
back to the dynamical calculation of ref. [55] . This calculation starts from an 
initial condition at r = 3.5 fm/c, which was given by a Parton Cascade Model 
(PCM) calculation for Au+Au collisions at y / Sj\T/v = 200 GeV. This initial 
condition is followed by a three dimensional Hubble flow. In addition, this 
dynamical calculation determined the thcrmodynamical constraints of entropy 
and local energy-momentum conservation for a sudden time-like deflagration 
from a supercooled quark-gluon plasma (QGP) phase to a pion gas phase. Such 
a transition may start at a constant proper-time, at some (position independent) 
value of the local temperature, that corresponds to the supercooled quark-gluon 
plasma (QGP) phase. Such sudden hadronization from a supercooled QGP 
phase can be realized due to the large characteristic nucleation times of hadronic 
bubbles inside a supercooled QGP, see ref. [21] for details. The strong three- 
dimensional expansion leads to negative pressures, mechanical instabilities and 
the resulting time-like deflagration may produce hadronic final states in the 
hatched area of figure 1 of ref. |25] ■ Half of this area corresponds to superheated 
hadron gas, with hadronic temperatures in the range of 1.0 - 1.4 T c , while half 
of available final states are hadron gas states below the critical temperature, 
with hadronic temperatures of 0.8 - 1.0 T c . Thus it is reasonable to assume, 
that in a realistic situation, the hadronic final states correspond to a range of 
temperatures even if the transition starts from the same temperature everywhere 
on a constant proper-time hypersurface. The variation can only be increased if 
the temperature of the prehadronic state happened to be inhomogeneous on a 
constant proper-time hypersurface. This is why we have assumed in the present 
paper, that the hadronic temperatures may be position dependent on a constant 
proper-time hypersurface. Instead of using a priori assumptions, we determine 
the parameter values of the hadronic local temperature distribution a posteriori, 
from recent data on Au+Au collisions at RHIC. 



3 Integration and saddle-point approximation 

The observables are calculated analytically from the Buda-Lund hydro model, 
using a saddle-point approximation in the integration. This approximation is 
exact both in the Gaussian and the non-relativistic limit, and if p v u v /T 3> 1 at 
the point of maximal emittivity. In this approximation, the emission function 
looks like: 

S(x, k)d±x = g gg^ffi) exp (-R-J(x - x.Y(x - x s y) d'x (17) 
27r d B(x.,p) + s a F 



G 



where 

Rrf = dnW-MSo))., (18) 
and x v stands here for (T,r x ,r y ,r z ). In the integration, a Jacobian r/t has to 
be introduced when changing the integration measure from d A x to drd 3 x. 
The position of the saddle-point can be calculated from the equation 

d fi (-]n(S ))(x s ,p) = 0. (19) 

Here we introduced 5*0, as the 'narrow' part of the emission function: 

s ^>'W^W q ' < 20 > 

In general, we get the following for the saddle-point in(r, ri x ,rj y ,rj z ) coordinates: 

Ts = TQ, (21) 

sinhr^ = ? P^ cosh ^ . (22) 



T J 1 7 ( a 

The system of equations (|22|) can be solved efficiently for the saddle-point po- 
sitions r] St i using a successive approximation. This method was implemented 
in our data analysis. For the widths of the distributions, we obtained exact 
analytic results, given in terms of rj Sy i that we determined from the successive 
approximations. We have used the full expressions, obtained simply with the 
help of derivations at the saddle-point, for the evaluation of the observables and 
for the comparison with the data. However, we summarize here only the leading 
order approximative result, due to reasons of clarity and transparency. 
In the simplest case, where all three rji s are small: 

Pi v 



Xi i+ * 



for i = x,y, z, (23) 
1 _ B{x s ,p) , , , . , (24) 



Ri,i B(x s ,p) + s q 

1 11 B(x s ,p) 




R 2 Ar| Ar2 B(x s ,p) + s Ar| ' 



(25) 



The above eqs. I|24l25ll imply, that the HBT radii are dominated by the smaller 
of the thermal and the geometrical length scales in all directions, as found in 
refs. Note that the geometrical scales stem from the density distribu- 

tion, governed by the fugacity term exp[ii(x)/T(x)}, while the thermal lengths 
originate from the local thermal momentum distribution exp[— p^u^ix) /T{x)\, 
and in the above limit they are 

1 m t d 2 



1 



R TA 




(26) 
(27) 
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4 The invariant momentum distribution 

The invariant momentum distribution can be calculated as 

Ni(p) = J d 4 xS(x, P ) = J d 4 xS c { P ,x). (28) 

The result is a simple expression: 

N ^ = ^ F ^ ^ y^P^ ' (29) 



where 



£ = p^(x s ), (30) 

F = { 2nf/^[detRl] 1/2 , (31) 

C 7 = 1 Zl (32) 
vA, i« 



Let us investigate the structure of this invariant momentum distribution, in 
particular, the exponent of the spectrum. Let us introduce 

b(x s , P ) =logB(x s ,p), (33) 

and evaluate this exponent in the limit, where the saddle-point coordinates are 
all small, 

6(*-,p) = 54r- + =4^ + 54r- + ^-^r-^ (34) 
^m t i* >x 2,m t l*. y 2m t ±*, z ±o 2m t i 

where the direction dependent slope parameters are 

T*, x = T + m t X 2 T °_ , (35) 
T + rata 1 

T*, y = To + m t Y 2 (36) 
T + m t a 2 

T*, z = T +m t Z 2 (37) 
T + m t a 2 

and 

m t = m t cosh(r) z . s - y). (38) 

In the limit when we neglect the possibility of a temperature inhomogeneity on 
the freeze-out hypersurface, a = 0, and using a non-relativistic approximation 
of Wit ~ m, we recover the recent result of ref. |27j for the mass dependence of 
the slope parameters of the single-particle spectra: 

T*, x = T +mX 2 , (39) 
T %y = Ta+mY 2 , (40) 
T,,, = To+mZ 2 . (41) 
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5 The elliptic flow 



Note, that b(x s ,p) is the only part of the IMD, that is explicitly angle dependent, 
so 



Ni(p) 



exp 





Pi ) 


2m t T* }X 












2m t T e ff 





Pi 



cos 



J ! 



where 



eff 



If J- 



(42) 



(43) 



So, we can easily extract the angular dependencies. Let us compute V2 by 
integrating on the angle: 

h{w) (44) 



where 



w 



l> 2 



Io(w) 



1 



4m t \T* >y 
Generally, we get from the definition 



1 

T 



Nl = 



d 2 n 



dpzPtdptdtp 27: dp z ptdpt 
the following equations: 



1 + 2 v n cos(nip) 



v 2n 



Injw) 

Io(w)' 



V2n+1 = 0. 



(45) 



(46) 



(47) 
(48) 



As first and the third flow coefficients vanish in this tilt angle $ has to 

be introduced to get results compatible with observations, as discussed in the 
subsequent parts. 

For large rapidities, |ry s — y\ becomes also large, and m t = wit cosh^ — y) 
diverges, hence w — > 0. Thus we find a natural mechanism for the decrease of 
V2 for increasing values of \y\, as in this limit, V2 — > li(0) / Iq(0) = 0. 



6 Elliptic flow for tilted ellipsoidal expansion 

Now, let us compute the elliptic flow for tilted, ellipsoidally expanding sources, 
too, because we can get a non- vanishing v\ and only this way, in case of i9 =/= 0, 
similarly to the non-relativistic case discussed in ref. ^ZT\. The observables are 
determined in the center of mass frame of the collision (CMS), where the r x 
axis points to the direction of the impact parameter and the r z axis points to 
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the direction of the beam. In this frame, the ellipsoidally expanding fireball, 
described in the previous sections, may be rotated. So let us assume, that we 
re-label all the x and p coordinates in the previous parts with the superscript ', 
e.g. x — > x' and p — > p' , to indicate that these calculations were performed 
in the system of ellipsoidal expansion (SEE), where the principal axis of the 
expanding ellipsoid coincide with the principal axis of SEE. In the following, we 
use the unprimed variables to denote quantities defined in the CMS, the frame 
of observation. 

We assume, that the initial conditions of the hydrodynamic evolution corre- 
spond to a rotated ellipsoid in CMS |27j . The tilt angle $ represents the rotation 
of the major (longitudinal) direction of expansion, r' z , from the direction of the 
beam, r z . Hence the event plane is the (r' x ,r' z ) plane, which is the same, as the 
{r x ,r z ) plane. The (zenith) angle between directions r z and r' z is the tilt angle 
while (azimuthal) angle ip is between the event plane and the direction of the 
transverse momentum pt . 

From the invariant momentum distribution, v m can be calculated as follows: 

2tt dn 

v m = / gggg co^mvQdy (49) 

We have made the coordinate transformation 

p' x = p x cosi? — p z sintf, (50) 

Py = Py, (51) 

p' z = p z cos$ +p x sin$, (52) 

and in addition: 

Px = pt cos ip, (53) 

p y = ptsinip, (54) 

and calculated the transverse momentum and the pseudorapidity dependence 
of 1)2, for a parameter set determined from fitting the axially symmetric ver- 
sion of the Buda-Lund hydro model to single particle pseudo-rapidity distribu- 
tion of BRAHMS [Hj and PHOBOS [§], the mid-rapidity transverse momentum 
spectra of identified particles as measured by PHENIX [23 EH! an( i the two- 
particle Bose-Einstein correlation functions or HBT radii as measured by the 
PHENIX [SI and STAR PS1 collaborations. 

We determined the harmonic moment of eq. I|49|l numerically, for the case of 
m = 2, but using the analytic expression of eq. I|29|l for the invariant momentum 
distribution, computing the coordinates of the saddle point with a successive 
approximation. The successive approximation means a loop here instead of 
solving the non-analytic saddle-point equations. We have chosen a loop long 
enough and have checked that an even longer loop will not modify the results. 
This was the same with the width of the integration-step. We integrated Ni(p) 
over p t , as the data were taken this way, too. Finally, we were able to describe 
112(77 = 0,Pt) and V2(r)) with the same set of parameters. 



10 



The results are summarized both in Fig. 1 and 2. We find that a small 
asymmetry in the expansion gives a natural description of the transverse mo- 
mentum dependence of V2- The parameters are taken from the results Buda- 
Lund hydro model fits to the two-particle Bose-Einstein correlation data (HBT 
radii) and the single particle spectra of Au + Au collisions at ^/snn = 130 GeV, 
ref. |15U26| . where the axially symmetric version of the model was utilized. Here 
we have introduced parameters that control the asymmetry of the expansion in 
the X and Y directions such a way that the angular averaged, effective source 
is unchanged. For example, we required that the effective temperature, T e g of 
eq. is unchanged. We see on Figs. 1 and 2 that this method was successful 
in reproducing the data on elliptic flow, with a small asymmetry between the 
two transverse expansion rates. 

The identified particle elliptic flow measurement of PHENIX used a method 
of determining the reaction plane from the particles at large rapidities, hence its 
results are not significantly affected by non-flow correlations, see ref. Fig. 1. 
illustrates the quality of agreement between our Buda-Lund model calculation 
and this PHENIX data set. From the parameter values corresponding to Fig. 
1, we calculate the value of the v-i{rj = 0) and find that this value is below the 
published PHOBOS data point at mid-rapidity by 0.02. Note, that in order to 
compute 1*2(77), one has to integrate Ni(r],p t , <fi) over p t first, and determine the 
elliptic flow from the p t integrated, rj and dependent spectra, as the PHOBOS 
data were taken this way, too. Because of this, there is no simple mathematical 
connection between Vzipti V — 0)j an d V2{v) as they do not stem from a common 
V2(j>t> v) function. 

The PHOBOS collaboration pointed out the possible existence of a non-flow 
contribution in their V2 data, see ref. j2j, as they did not utilize the fourth order 
cumulant method to determine v%. We attribute the 0.02 difference between 
the present Buda-Lund model calculation and the PHOBOS data point at mid- 
rapidity to such a non-flow contribution |34|. The magnitude of the non-flow 
contribution has been explicitly studied (but in a different acceptance, at mid- 
rapidity) by the STAR collaboration. STAR found that its value is of the order 
of 0.01 for mid-rapidity minimum bias data in the STAR acceptance, ref. p]. 

The good description of the dn/drj distribution by the Buda-Lund hydro 
model |151 l2l)| is well reflected in the good description of shape of the pseudo- 
rapidity dependence of the elliptic flow. Thus the finiteness of the expanding 
fireball in the longitudinal direction and the scaling three dimensional expansion 
is found to be responsible for the experimentally observed violations of the boost 
invariance of both the rapidity distribution and that of the collective flow Vi- 

Furthermore, the parameter values corresponding to Figs. 1 and 2 indicate a 
high, T > T c = 170 MeV central temperature, with a cold surface temperature 
of T s sw 105 MeV. The success of this description suggests that a small fraction 
of pions may be escaping from the fireball from a superheated hadron gas, which 
can be considered as an indication, that part of the source of Au + Au collisions 
at RHIC may be a deconfined matter with T > T c . 

Let us determine the size of the volume that is above the critical temperature. 
Within this picture, one can find the critical value of s — s c from the relation 
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that r /(l + as c ) = T c . Using T = 210 MeV, T c = 170 MeV, and a = 1 wc 
find s c = 0.235. The surface of the ellipsoid with T > T c is given by 

t--2 ^2 2 

W+# + #= L (55) 

The principal axes of the "critical" ellipsoid are given by X c — Xf^fs^ ~ 4.2 
fm, Yc = Yfy/s^ ~ 5.1 fm, Z c = Zf^fs^ ~ 8.5 fm, hence the volume of the 
ellipsoid with T > T c is V c = 4fA c y c Z c w 753 fm 3 . 

Note, however, that the characteristic average or surface temperature of the 
fireball within this model is T s = To/(l + a) « 105 MeV. So the picture is similar 
to a snow-ball which has a melted core inside. 

Our study shows that this picture is consistent with the pseudorapidity and 
transverse mass dependence of t>2 at RHIC in the soft pt < 2 GeV domain, 
however, it is not yet a direct proof of the existence of a new phase. Among 
others, we have to determine precisely the errors on the best fit parameters and 
to determine the confidence levels of the fits, which will be a subject of further 
research. 

In discussing the significance of the results, it is useful to compare it to other 
calculations to find similarities and differences as well as to map out directions 
for possible further research. For us the key point is not the good agreement 
between the model and the data, but the analytic insight and the functional 
relationships between the model parameters and the observables. We checked 
the model against the data only to demonstrate that we are on a good track to 
understand the rapidity dependence of the elliptic and higher order flows, but 
the fine-tuning of the model parameters is a subject of further investigations. 

Our most important result seems to be eqs. (|38I45I47I48I) . that explain ana- 
lytically why all higher order flows vanish at very forward or backward rapidities. 
Due to the finite longitudinal size of the source, the point of maximal emittivity 
moves to the summit of the expanding ellipsoid as the rapidity is increased to 
high values. Due to the Hubble flow, the transverse momentum distribution has 
negligible transverse flow contributions at this point, the local temperature plays 
the dominant role. However, the local temperature contributes equally in both 
transverse directions, hence all second and higher order flows vanish at very 
forward rapidities. Similar observations hold in the very backward direction, 
due to symmetry reasons. 

When comparing to earlier calculations, we observe that the two key features, 
the finite longitudinal size and the three dimensional Hubble flow were not 
present simultaneously in other works as far as we know. Also, the temperature 
variations, the cooling on the surface were not considered by other attempts to 
understand elliptic flow. For example, Hirano and Tsuda considered a three- 
dimensional numerical solution of relativistic hydrodynamics in ref. |35j . Their 
model is not too far from the considerations presented here, they have a finite 
longitudinal extension of the source and a well developed transverse flow at 
mid-rapidity. Their Fig. 9 indicates that they obtained a vanishing elliptic flow 
at very forward and backward rapidities. However, they utilized a Bjorken type 
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of initial condition, and the concept of a constant freeze-out temperature. As a 
result, their elliptic flow is too flat, too Bjorken-like near mid-rapidity. Hirano 
studied in ref. the effects of short lived resonance decays on the rapidity 
dependence of the elliptic flow at SPS energies. His result is that resonance 
decays yield a negative non-flow contribution, ranging from -0.15 at mid-rapidity 
to about -0.05 at forward and backward rapidities. We did not explore the 
consequences of such an effect here, however, its magnitude is about the size of 
the error bars on the PHOBOS data. This is one of the interesting directions 
that can be explored in further studies, and the importance of this effect will 
increase at RHIC as more and more high statistics, precision measurements will 
be available on the elliptic flow. 

7 Summary and conclusions 

We have generalized the Buda-Lund hydro model to the case of ellipsoidally 
symmetric expanding fireballs. We kept the parameters determined from fits to 
the single particle spectra and the two-particle Bose-Einstein correlation func- 
tions (HBT radii) and interpreted them as angular averages for the 
direction of the reaction plane. Then we found that a small splitting between 
the expansion rates parallel and transverse to the direction of the impact pa- 
rameter, as well as a small zenith tilt of the particle emitting source is sufficient 
to describe simultaneously the transverse momentum dependence of the collec- 
tive flow of identified particles £Q at RHIC. If a constant non-flow parameter of 
0.02 is added to the calculated values, we also describe well the pseudorapidity 
dependence of the collective flow 0E]. 

The results confirm the indication for quark deconfincmcnt at RHIC found 
in refs. |15U26| . based on the observation, that some of the particles are emitted 
from a region with higher than the critical temperature, T > T c = 170 MeV. 
We estimated that the size of this volume is about 1/8-th of the total volume 
measured on the t = tq main freeze-out hypersurface, totaling of about 753 
fm 3 . However, the analysis indicates that the average or surface temperature is 
rather cold, T s w 105 MeV, so approximately 7/8 of the particles are emitted 
from a rather cold hadron gas. 
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Figure 1: Buda-Lund model calculation and the ^(ptj V = 0) data 
Ellipsoidally symmetric Buda-Lund calculation compared to the PHENIX V2(pt) data 
of identified particles Q. The parameter set is: To = 210 MeV, X = 0.57, Y = 0.45, 
Z = 2.4, a = 1, T = 7 fm/c, ■& = 0.09, X/ = 8.6 fm, Y; = 10.5 fm, Z f = 17.5 fm, 
{Jiq.ty — 70 MeV, ho,k = 210 MeV and /io,p = 315 MeV, and the masses are taken as 
their physical value. 
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Figure 2: Buda-Lund model calculation and the V2(tj) data 
This image compares the ellipsoidally symmetric Buda-Lund model to the 130 GeV 
Au+Au and 200 GeV Au+Au v 2 (r]) data of PHOBOS 0111 Here we used the same 
parameter set as at fig. with pion mass and chemical potential, and added a constant 
non-flow parameter of 0.02 to the calculated values of v 2 . 
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Parameter 


Value 


To 


210 MeV 


X 


0.57 


Y 


0.45 


Z 


2.4 


a 


1 


TQ 


7 fm/c 


At 


fm/c 




0.09 


Xf 


8.6 fin 


Yf 


10.5 fm 


Zf 


17.5 fm 


A*0,TT 


70 MeV 




210 MeV 


A*0,p 


315 MeV 



Table 1: The table shows the parameter set used to describe the data. Note, 
that a 2 = T ° 7 7 Ts , see equation JTBJ. 
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